Solitons in cavity-QED arrays containing interacting qubits 
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We reveal the existence of polariton soliton solutions in the array of weakly coupled optical cavities, 
each containing an ensemble of interacting qubits. An effective complex Ginzburg-Landau equation 
is derived in the continuum limit taking into account the effects of cavity field dissipation and qubit 
dephasing. We have shown that an enhancement of the induced nonlinearity can be achieved by two 
order of the magnitude with a negative interaction strength which implies a large negative qubit-field 
detuning as well. Bright solitons are found to be supported under perturbations only in the upper 
(optical) branch of polaritons, for which the corresponding group velocity is controlled by tuning 
the interacting strength. With the help of perturbation theory for solitons, we also demonstrate 
that the group velocity of these polariton solitons is suppressed by the diffusion process. 



I. INTRODUCTION 

To create new devices for quantum information stor- 
age and processing, macroscopically coherent coupled 
matter-light states are at the heart of researches, see e.g. 
[l|. Cavity quantum electrodynamics (cavity-QED), the 
study of coherent interactions between matters and quan- 
tized cavity electromagnetic fields, have provided a versa- 
tile and controllable platform to describe many interest- 
ing phenomena in this field; they are micro-masers, Pur- 
cell effect, squeezed state generation, atom trapping, and 
quantum state transfer, see e.g. The simplest light- 

matter system, one cavity mode interacting with a sin- 
gle two-level atom, is described by the Jaynes-Cummings 
model [J]. As the number of two- level atoms increases, 
collective effects, due to the field mediated interactions of 
atoms among themselves, can be described by the Dicke 
Hamiltonian and give rise to intriguing many-body phe- 
nomena [5|. Apart from a single cavity configuration, 
cavity-QED arrays composed by engineered optical cav- 
ities modes, few-level atoms and laser light modes, may 
serve as a many-body system for light [6| . Novel quantum 
phase transitions of light, from Mott-insulator, glassy, 
to super-solid states have been demonstrated to analyze 
critical quantum phenomena in conventional condensed 
matter systems by manipulating the interaction between 
photons and atoms 

From practical point of view, two important prerequi- 
sites seem to be useful to provide quantum information 
processing and transmission. In particular, interactions 
are essential to generate necessary quantum correlations 
and quantum entanglement [l|. Second, weakly coupled 
cavity-QED arrays and lattices represent suitable candi- 
dates for quantum information transport by using pho- 
tonic tunneling effect In particular, it is worth notic- 
ing possible application of diamond NV centers for cou- 
pled cavity-QED arrays [12]. The qubits based on NV 
centers are potentially attractive because they are posing 
essentially longer dephasing time (up to few milliseconds) 
even at high enough temperatures [l3j. 

Another comprehensive example on this way is the res- 



onant excitation of a single quantum dot (QD) strongly 
coupled with a photonic crystal nanocavity experimen- 
tally demonstrated in |14l - [la |. Moreover, measurement of 
time-dependent second order autocorrelation function in 
such experiments enables to distinguish photon blockade 
regime that is connected with anharmonic energy-level 
spacing and is, obviously, relevant to nonlinear effects in 
QD-light interaction in small volume cavities. Noticing 
that the detunings between frequency of coherent probe 
beam (or frequency of resonant exciton state), and the 
bare-cavity frequency represent vital parameters in such 
experiments, cf. Til 16j. 



In the future such experiments permit to create on-chip 
coupled cavity-QED arrays with various architecture for 
few photon light field processing and transport integrated 
with telecom lines |17| . It is known that quantum optical 
solitons are a natural candidate for quantum optical in- 
formation prosessing. Studies of their transmission, for- 
mation, and transformation in coupled cavity-QED ar- 
rays represent important steps for achieving this aim. 
Polariton solitons have some advantages in this case. In 
particular, as demonstrated in [18J, polaritonic nonlinear- 
ity can be high enough in comparison with pure optical 
nonlinearities achieved, e.g. with pure optical cavity soli- 
tons in VCSELs [l!| , and can permit to achieve a soliton 
regime for essentially smaller particle (polariton) num- 
ber. 

More precisely, here we deal with the problem of po- 
lariton soliton formation occurring at quantum matter- 
filed interface in semiconductor quantum wells (QWs) 
embedded in Fabry-Perot microcavities. Nowadays po- 
lariton solitons and relevant superfluid behavior of non- 
equilibrium exciton-polaritons with narrow-band (GaAs) 
semiconductor structures are observed in a few labs, [l8|, 
[2fj| . The main physical features of such solitons are 
connected with the balance between dissipation ef- 
fects and external pumping that occurs for parametri- 
cal processes allowing to create non-equilibrium exciton- 
polariton Bose-Einstein condensate (BEC) at lower po- 
lariton branch under the temperatures of few kelvins. 
The solitons in experiments are excited within picosecond 
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time scale and localized in micrometer-scale sizes. 

Recently, in [2l[ we proposed a 2D polaritonic crystal 
model that represents cluster material exhibiting high 
nonlinear properties due to the small but macroscopical 
number of two-level atoms strongly coupled with opti- 
cal field in the cavity lattice. Basing on the Holstein- 
Primakoff approach, we have shown that such nonlinear- 
ity can basically be created by means of two physical 
processes: that is a saturation effects occurring due to 
localization of macroscopically small number of atoms at 
each cavity, and atom-atom interaction on its own. Ex- 
perimentally such a system can be designed by using 2D 
photonic crystal host with defect cavities [22j . Noticing 
that exciton-photon saturation effects can also be strong 
enough and accompanied by Coulomb repulsion processes 
for coupled cavity QED arrays based on the appropriate 
number of multiple QWs structures, cf. (23[. Coherent 
matter-field interaction can be achieved in this case at 
room temperatures by using wide-band (ZnO) semicon- 
ductor microstructures, cf. [24] . 

In the paper we consider the problem of soliton forma- 
tion in ID weakly coupled cavity arrays containing two- 
level systems (atoms, QDs, etc) that interact with sin- 
gle mode optical cavity fields. Instead of non-interacting 
(free) particles, several systems support interacting com- 
ponents which are referred to the qubit for simplicity. For 
the ensemble of interacting qubits, the Lipkin-Meshkov- 
Glick (LMG) model is an example of the system produc- 
ing maximal pairwise entanglement under the phase tran- 



sition of its ground state [25[. Recently, it has been pro- 



posed that the LMG model, originated in nuclear physics, 
may describe Josephson effect or two-mode Bose-Einstein 
condensate (BEC) (2^]. In particular, relevant nonlin- 
ear model which is called in [27j as " canonical Josephson 
Hamiltonian" evokes great interest in the framework of 
efficient spin squeezing [28] and resent experiments with 
high precision measurements [29| . Besides, the schemes 
for realizing a dissipative LMG model in optical cavity- 
QED 30( and in circuit-QED 3l| have been discussed. 
By introducing the interaction with a quantized optical 
cavity mode, the ground state phase transition [32[ , max- 
imal shared bipartite concurrence of the ensemb le 1331 . 
and dynamics from Rabi to Josephson oscillations 34j,|35[ 
for the extended Dicke-LMG (DLMG) model have been 
studied. 

In this work, we extend DLMG model to the array of 
quantized field cavity modes and study nonlinear dynam- 
ics in such a qubit-cavity-QED array system. In the con- 
tinuum limit, we use a multiple-scale envelope function 
(MSEF) method (3(| H?} to construct soliton solutions 
for considering nonlinear interactions among qubits as 
well as dissipation effects for them. In particular, we de- 
rive and examine a complex Ginzburg-Landau equation 
(GLE) that supports "slow" soliton propagations in the 
presence of polariton formation. Being a step in the stud- 
ies of collective properties of light with interacting media, 
our results in the proposed qubit-cavity-QED system are 
shown to be connected with the formation of coupled 
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FIG. 1. (Color online) Schematic for our proposed ID qubit- 
cavity QED arrays, in which each cavity contains an ensemble 
of qubits with the interaction Hamiltonian among them de- 
scribed by the LMG model, -Hlmg oc Jf . 



matter-field states, upper (UB) and lower (LB) branch 
polaritons, which would be natural carriers for quantum 
information processing at matter-field interface. 



II. MODEL OF INTERACTING 
QUBIT-CAVITY-QED ARRAYS 

As illustrated in Fig. Q] we consider a one-dimensional 
(ID) array of optical cavities, each containing an en- 
semble of Ni of interacting qubits. The interaction in- 
side each single cavity is described by Dicke-LMG model 
Hamiltonian, Hf hMG , and is represented as, cf. [27r - t30lj 
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Subscript i labels the i-th cavity to be used afterwards. 
Detuning frequency Aj has been defined as the difference 
between qubit transition frequency Wj and quantized field 
frequency uif, i.e., Aj = CJj — uif, orbital angular mo- 
mentum representation, Jj^ with j — x, y, z, ±, has been 
chosen to describe the qubit ensemble, while the field 
is described by its creation (annihilation) operators, aj 
(a,); the field-ensemble and the inter-qubit coupling are 
represented by the constants <?j and rji. 

Hamiltonian in Eq. ((!} is equivalent to the LMG 
model in the limit g,; = 0, and hereby is called Dicke- 
LMG Hamiltonian, which can be obtained from the 
Gross-Pitaevskii equation describing a two-species BEC 
(see e. g. 27]) interacting with a quantum field [3^ |. 
Experimental realizations providing an assorted range 
of tunable parameters for the DLMG model may in- 
clude a two-hyperfine-structure-defined-modes BEC cou- 
pled with a quantum cavity field mode through one mi- 
crowave photon process; e.g. trapped hyperfine ground 
states of a Sodium BEC inside a microwave cavity [33| . 
In the area of solid state physics, superconducting qubit 
devices seem to be promissing candidates to realize qubit- 
cavity-QED being under discussion, see e. g. (38[- In par- 
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ticular, arrays interacting superconducting qubits cou- 
pled with the quantum field mode of a coplanar waveg- 
uide resonator may be provided the ensemble sizes are 
small [H. 

For the configuration of an array of cavities, the qubits 
are assumed to be confined inside the cavities, while pho- 
tons can travel from one cavity to another. We also 
assume that the traveling process for photons is domi- 
nated by the hopping between adjacent sites. The value 
of hopping matrix elements is a function of the distance 
between cavities. If the distance between each cavity is 
sufficiently small, a cavity field will only hop to the near- 
est neighbors. Then the total Hamiltonian of interaction 
for this qubit-cavity-QED array is given as (cf. 0), 
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where a is a photon-hopping strength from cavity site i 
to i + 1, and M counts the number of cavities. In the 
present work, we assume that only one species of qubits 
is considered and that the number of ensemble qubits 
is the same in each cavity, i.e., Tji = 77, u>{ = lo, and 
Ni = N. For simplicity, the parameters are taken as 
real and homogeneous so that the detuning frequency, 
hopping and coupling strengths be all identical in each 
cavity site, i.e., Aj = A, gi = g, and cti = a, respectively. 
Physically, this approach is valid if the separation of two 
adjoint cavities is macroscopically large in comparison 
with the atomic de Broglie wavelength and overlapping of 
atomic wave functions from neighbor cavities is negligibly 
small. 

If most of qubits stay in the ground state, we can 
apply Holstein-Primakoff transformation [39| by map- 
ping qubit excitation operator Jjj on the Schwinger 
representation [4(| for a two-level system, i.e., J+j — 



Nb! ( 1 



b\bi/N) and J z 
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with the commutation relation [bi,bj] = Sij. Practi- 
cally, the number of qubits N is large enough and one 
can expand the orbital angular momentum representa- 
tion for the qubit ensemble up to the order of 1/N, i.e., 
[1 - {blk/N)} 1 / 2 rj 1 - b\bi/2N. We are working within 
the so-called low-excitation-density limit when the qubits 
mostly populate their ground levels. 

Since a photonic field is an order parameter for the sys- 
tem described by Hamiltonians ([T]), ([2]) it is non-zero for 
superfluid phase of coupled matter-field state in the pres- 
ence of photon tunneling between neighbor cavities, see 
e.g. At the same time low-excitation-density limit 

implies that the average photon number at each cavity 
should be essentially smaller than the average number of 



a] a, 



qubits at the excited state, that is 

cf. mi. 

Now, since number N is assumed to be large enough, 
we are able to neglect higher-order terms in the expansion 
of J + j operator. Then, we can transform the Hamilto- 



nian in Eq. ([2]) into 
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Based on this Hamiltonian, Heisenberg's equations of 
motion for the involved two bosonic operators of cavity 
modes on and qubit excitation bi are 
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r/bi + -j-fblbih. 

(5) 



Here, we have neglected quantum fluctuations due to 
a large number of qubits considered. 

To have mean-field solutions, we replace the pair of 
operators by the corresponding expectation val- 

ues and approximate this array configuration 

to a continuous model, i.e., Vi+i + V^-i ~ 2tp(x,t) + 
a ^sip{x, t), with the distance between two adjoint op- 
tical cavities denoted by d. By considering the conserva- 
tion of the total photon number and qubit excitations, we 
can re-normalize the variables with respect to the num- 
ber, iVpoi, defined as iV po i = |?/;| 2 + |/3| 2 . The conserved 
quantity ./Vp i is also known as the number of polaritons 
at each cavity (42|. By replacing the cavity field and 
qubit excitation by -0 — > ip/y/N po i and (3 — > fi/-JN po \, 
respectively, now we have a normalization condition for 
the cavity field and qubit excitation, \ip\ 2 + \/3\ 2 = 1. 

Thus, in this continuum limit, the equations of motion 
in Eqs. (14151) become, 



iftV = (A — 2a - i7c )^ - ad 2 d xx ^ + g/3 - £/ sat |/3| 2 /3, 



id t /3 = -(iT d + ri)l3 + gi> 

- U sat [p 2 ^ + 20|/3| 2 ] + C/ int |/3| 2 /3, 



(6) 
(7) 



where U int = 2r]n po \, U sat 



g" P oi . 



ipoi = N po i/N is a 



polariton number density. In Eqs. (|6I7[) we have a phe- 
nomenologically introduced decay rate for photonic field 
7 C that characterizes the leakage of photons in the cavity 
and dephasing rate for the qubit system. Physically, 
parameter J7 sa t characterizes a nonlinear saturation effect 
under the qubit-light interaction. Parameter U{ n t is re- 
sponsible for qubit-qubit interactions within the LMG 
model. Let us briefly discuss the applicability of our 
model (Eqs. (|6I7I0 for different physical systems repre- 
senting qubits. 

Example 1: cavity QED with atomic qubits. Here we 
discuss cavity QED array with two-level atoms as a qubit 
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system. To be more specific, we consider ultracold two- 
level rubidium atoms with resonance frequency u; Q f,/27r = 
382THz that corresponds to mean weighted rubidium D- 
lines (42[. The atomic polarization dephasing rate and 
the minimal value of each cavity field decay rate can be 
taken as several tens of megahertz 's that corresponds to 
cavity quality factor Q = ui a i,/2j c ~ 10 6 . The strength of 
interaction of a single atom with a quantum optical field 
is taken as go/2ir — 89.5MHz at each cavity with the 
effective volume of atom- field interaction V — 5000/im 3 . 
To achieve a strong atom-field coupling regime, see Eq. 
(|15p below, one can propose a macroscopically large num- 
ber of atoms at each cavity, say N = 5 x 10 5 . This num- 
ber is relevant to the density p = 10 14 cm~ 3 of ultracold 
atoms that implies a collective atom-filed coupling pa- 
rameter as g = goVN — 2tt x 63.2GHz at each cavity. 
At the same time the parameter that describes atom- 
atom interaction n = 47rha sc p/m in the Born approxi- 
mation can be estimated as rj — 2n x 0.73kHz, where a sc 
is atomic scattering length that we choose as a sc = 5nm, 
cf. Thus, for atomic QED array U int /U sat ~ 1CT 8 

and we can neglect atom-atom interaction processes in 
Eqs. (07]). 



III. DISPERSION RELATIONS AND GROUP 
VELOCITIES 

To construct solitons in this cavity-QED arrays with 
interacting qubits, we seek for wave-packet solution s by 
the multiple-scale envelope function method [36], l37j |. 
With the introduction of different length and time scales, 
i.e., x m = \ m x (A < l,m = 0, 1,2...) and t m = \ m t 
(A < l,m = 0,1,2...), we can expand photon field, ip, 
and qubit excitation, /3, as 

i\> = XiP (1) + \ 2 4j (2) + A 3 ^ 3 ' + . . . 
P = + \ 2 f3 (2) + \ 3 P {3) + ... 

By substituting these expansions in Eqs. (JB][7]), we can 
gather all terms that are proportional to the first order 
of A and find, 

id i0 V (1) = (A - 2a - i 7c )^ (1) - ad 2 d 2 ip {1) + <?/3 (1) , 

(8) 

id t0 P {l) =-(iT d + ri)pW +g^ l) . 

(9) 



Example 2: cavity QED with excitonic qubits in QWs. 
A 2D version of Eqs. (|6I7|) describes properties of exciton- 
polaritons that can be created in a microcavity with 
Bragg mirrors under strong coupling between excitons 
confined in QWs and a cavity photonic mode, see 
e -g- [HHHl- It has been shown recently that the strong 
coupling condition can be achieved in the system of pho- 
tonic crystal polaritons (45|, which can be a possible 
platform to implement qubits being proposed in a cavity- 
QED-array. The typical values of Rabi splitting fre- 
quency for QWs based on GaAs or CdTe is about few 
of meVs. The saturation term U sa t in Eqs. (|6l7p for 
QWs structure can be represented as U sa t — g/n sa tS, 
where n sat is exciton saturation density, S is a quanti- 
zation area. The typical values for n sat and U sa t are 
1.4 x 10 n cm~ 2 and 2ir x 0.65MHz, respectively, taken 
for exciton Bohr radius lOnm and excitation laser spot 
S ~ 400/im 2 . The exciton-exciton interaction term in 
([7| is evoked by Coulomb interaction and about two or- 
der times larger than U sa t, cf. [44|]. Thus, we can neglect 
exciton-photon saturation effects for excitonic qubits. 
Typical experimentally accessible values for relaxation 
processes in semiconductor structures can be approached 
as Td — 2n x 12.1GHz, j c ~ 2n x 50GHz that implies 
exciton-polariton lifetime within tens of picoseconds. 



This first-order expansion for the qubit excitation j3^ 
in Eq. §§§ supports a plane wave solution, from which 
we can find a corresponding dispersion relation by using 
the solution in the form of = E (1) e^o-^o) an d 



rier frequencies are 
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1/21 



(10) 



where k is a wave number. Physically, Eq. (|10[) repro- 
duces a familiar result for frequencies of two branches of 
polaritons [46j|, which are denoted as UB (w+) and LB 
(cj_), respectively. 

It is important to note that field (ip^) and qubit 
excitation (fi^ 1 ') variables become exact solutions of 
Eqs. (|6l7p neglecting nonlinear effects, i.e., setting U sat = 
Uint — 0. Hence, it is easy to understand that MSEF 
method developed here is valid if dispersion character- 
istics of polaritons, see Eq. (fT0|) . cannot be essentially 
modified taking into account nonlinear effects. Strictly 
speaking, the condition 



The set of coupled nonlinear equations (I6I7[) is the 
starting point of our work, and soliton solutions both 
for the wave-packet envelope of cavity field ^ (order pa- 
rameter) and for qubit excitations /3 are considered ana- 
lytically below. 



U sat ,U int < 2g (11) 

should be fulfilled for coupled qubit-light systems de- 
scribed by Eqs. (|6I7[) . In practice this condition implies 
achieving low-excitation density limit n po \ <ti 1 discussed 
above. 
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Next, for the second order of multiple-scales, i.e., A 2 , 
we can have a linear wave equation for a wave-packet 
envelope 



(d tl +v ± d Xl )E^ =0, 



(12) 



from which one can find the corresponding group veloci- 
ties, 



v± = d k uj ± = 



2akd 2 Q% 
^±+9 2 



(13) 



that are defined for UB (u + ) and LB (u_) polariton wave- 
packets. Below we consider variable in connection 
with Eq. (I12[) in the moving frames £± = Xi—v±t%. In Eq. 
([TBf we have also introduced characteristic frequencies 



n± = - 

* 2 



S-i<y± ({8- iif +Ag 2 \ 



1/2 



(14) 



where 7 = 7 C — Yd is an effective qubit-field decay rate for 
the system and 5 = A +77— 2a (l — d 2 k 2 /2) is a total mo- 
mentum dependent qubit-light detuning. Now, the qubit 
interaction parameter 77 introduces an additional phase 
shift that can be used to tune the relative dispersion re- 
lations of polaritons. 

It is noted that the group velocities v± defined in Eq. 
(|13[) are complex numbers in general. It is known that 
this leads to the additional deformation of a pulse enve- 
lope propagated in the medium [43|. The above MSEF 
method works only for small decay rates if a strong cou- 
pling qubit-light condition is valid, that is 



7c T d < 2.9. 



(15) 



In particular, in this limiting case the imaginary part 
of frequencies Im(a;±), which characterizes the decay of 
polaritonic systems, is smaller than the real one Re(w±), 
i.e., Im(w±) <C Re(w±). If the condition in Eq. (|T5|) is 
not met, the group velocity can be correctly determined 
in terms of the energy flow in the medium only [48] . How- 
ever, this is not the case of our present consideration. 

With the introduction of the qubit interaction, in Fig. 
2(a), we show the dependence of these two branches 
10 ± on the wave number k under the resonance condi- 
tion A = 0, but for different qubit interaction strengths, 
77 = —1, 0, and 1, respectively. Equation (fTUf is valid for 
a small quasi-momentum k, such as kd -C 1. In particu- 
lar, if we approximate 1 — d 2 k 2 /2 w cos(fcd) in Eq. (ITUl) . 
then, by setting 77 = and 7 C = Td = 0, we can recognize 
true dispersion relation taken at first Brillouin zone for 
periodic structure, see Fig. 2 and cf. [2l[. For a small 
quasi-momentum k we can have parabolic dispersion re- 
lation that is inherent in LB polaritons occurring in semi- 
conductor microcavities [Ullloj]. From Fig. 2(a), one can 
see that a negative value of 77, accounting a repulsive in- 
teraction among the qubits, moves the two branches u>± 
upward to a higher value; while a positive one moves the 
dispersion relation downward. 
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FIG. 2. (Color online) (a) Two branches in the dispersion 
relation, labeled as upper lu+ and lower ui- branches, are 
shown as a function of the wave number k at the resonance 
condition A = 0, but for different qubit interaction strengths, 
77 = —1, 0, and 1, respectively, (b) The carrier frequencies w± 
and (c) the group velocity v± for the wave packet are shown 
as a function of interacting strength 77 under the condition 
k — 0.1, 7 = 0, but with different frequency detunings, i.e., 
A = —2, 0, and +2, respectively. Other parameters used are 
a = d = 1, g — 10, and n po i = 0.01. 
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By collecting terms of ij — A in Eq. (JTOj) under the con- 
dition (|15p , we can make the linearization approximation 
for the dispersion relation, i.e., 



a I cos(fcd) 



V 



2a 



±\9 



(16) 



where one can view this qubit interaction as a 
linear frequency detunging under the conditions of 

2Aacos(kd) < a 2 (cos(fcd) + <| g | 2 . In this 

case, we obtain the dispersion relation as a function of 
r] — A only, i.e., the nonlinear interaction strength (77) is 
such a linear term as photon-atoms detuning factor (A). 
Fig. 2(a) shows such a shifting of dependencies. 

By fixing a wave number k = 0.1, in Fig. 2(b), the 
dependences of these two frequency branches are mod- 
ified strongly with respect to the change of interaction 
strength 77. As the interaction energy among the qubits 
changes from a repulsive to an attractive one, both fre- 
quency w + and w_ decrease as the value of the qubit 
interaction strength increases. But when the above lin- 
ear condition of r\ — A is not valid, these modified fre- 
quency curves for two branches demonstrate a totally 
different tendency to the frequency detuning, A. For a 
large enough negative (positive) value of the interaction 
strength, 77 -C —10 (77 ^> 10), the upper (lower) branches 
for different detunings merge into a single curve. That is, 
in this condition we have a frequency curve which is inde- 
pendent on the detunings. Now a nonlinear interaction 
term 77 is dominant. 

An example of the group velocity v± for two branches 
is shown in Fig. 2(c) as a function of qubit interaction 
strength 77, for k = 0.1 under the condition (fl~5|) but 
with three different frequency detunings, i.e., A = —2, 0, 
and +2, respectively. It is noted that the group velocity 
in this work is normalized to ad, instead of the speed 
of light. Behaviors of the group velocity of light in the 
cavity array under discussion can be easily understood 
by using a physical picture of polariton propagation - 
see Sec. IV. Here we note, that the difference between 
these group velocities in the two branches becomes larger 
while the absolute value of the qubit interaction strength 
77 increases. A "slow" propagation of the wave-packets 
can be achieved just by changing the qubit interaction 
parameter 77 for a negative (positive) value for the UB 
(LB) solution. 



IV. POLARITON SOLITON SOLUTIONS 

Now, we do a multiple-scale expansion up to the third 
order by collecting terms with A 3 , and reach a complex 
GLE El, i.e., 



"at 



d 2 ^ 
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(17) 



where variables * = A£ (1) , X± = £±/A, and t 2 = X 2 t 
have been introduced. The coefficients D± and C± , ap- 



pearing in the second and third terms of Eq. ([17]) respec- 
tively, have the forms, 



D+ = 



ad 2 n* 



v 2 ± g 2 



(n 2 ± 



c± 



ooi9 4 [(30t + fi£) - 477] 



2 0, 



(18) 
(19) 



Coefficients defined in Eqs. (|18|19|) are complex and can 



be evaluated as D 



± 



D 



(i) 



c 



(1) 



(2) 



Real part D ± of D±-coefficient describes diffraction ef- 
fects occurring with the wave-packet; while the imaginary 
part D± characterizes diffusion processes. Parameter 

is responsible for a Kerr-like nonlinearity that oc- 
curs due to polariton-polariton scattering. At the same 

time an imaginary part, that is C± , is relevant to non- 
linear absorption effects. 

By taking into account Eq. (ITU)) , it is helpful to intro- 
duce convenient UB (m+) and LB (m_) polariton masses 



m± 



H [dlu± 



m ph 



g 



fe=0J 



n 2 ± 



(20) 



where m v h — h/2ad 2 is an effective photon mass in the 
cavity [411 ]. Neglecting the kinetic energy of polaritons 
and supposing that the solutions are taken at the bottom 
of dispersion curve, one can approximate 



2m± 



C± 



277 D0 lOy 3 , 

± Y 1 (21) 



Y, — JL 



1± 



where we have introduced Hopfield coefficients such as 
, s and for the sake of simplicity 

y/S 2 +ig 2 J 

we suppose that g > 0. The latter parameters determine 
the contributions from photons and qubit excitations to 
polaritons. In particular, we can express field variable 
\1/ that is a flight qubit through polariton qubit states as 
following 



UB 



LB, 



(22) 



where ^ub and Elb are new variables characterizing UB 
and LB polaritons respectively. They are defined in a 
convenient way by using Bogoliubov transformation [2l| . 

In particular, for a positive and large frequency de- 
tuning 8, that is \S\ g and S > 0, we have Y + s» 1 
and y_ w t|t, respectively, which corresponds to photon- 
like UB polaritons (*E,ub — ^) with mass m+ ss m p h 
and group velocity v + = hk/m + , see Eqs. (|20H22[) and 
Fig. 2(c). Remarkably, at the same time the group veloc- 



ity of LB polaritons with mass m_ 



approaches 



V— = hk/m,— <C c, where c is speed of light in vacuum. 
Now, we have a "slow" (matter-like) soliton formation in 
the cavity array. For \5\ 3> g but S < 0, a physical picture 
becomes opposite and a propagating optical pulse rele- 
vant to LB polaritons (E.lb — with mass to_ w m P h 
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FIG. 3. (Color online) Nonlinear coefficient C± of the reduced 
NLSE in Eq. (|17p is shown for two frequency branches uj± , as 
a function of the qubit interaction strength g for different 
frequency detunings A. Inset shows a blowup close to g = 0. 
The parameters used are k — 0.1, a — d — 1, n po \ = 0.01 and 
g = 10, respectively. 



and group velocity i>_ = hk/m— 3> w+. In the presence 
of qubit-light resonance condition, 6 = 0, UB and LB 
polaritons equally contribute to the qubit state, that is 
Y± = 1/ The masses of polaritons are equal to each 
other, i.e., m po i = m± = 2m p h- Flight qubit propagates 
in the medium with the velocity v± = hk/m po \. 

At thermal equilibrium (or quasi-equilibrium), the 
lower polariton branch is macroscopically occupied. In 
this case we are interested in the properties of LB po- 
laritons only. Taking into account the dependence of ^ 
on LB polariton variable Slb in Eq. (|22[) . it is easy to 
rewrite the GLE in Eq. (IT7|) for LB polariton variable as 



LB 



h d 2 E 



LB 



Of 



2m_ dXl 



C p \^lb\ 2 Slb 



0, (23) 



where parameter C p ~ 2n po \gY?/Y_ establishes two- 
body polariton-polariton interaction strength. Thus, 
Eq. ([23| derived above by MSEF method remains in 
complete agreement with the results obtained for LB po- 
laritons by using Hamiltonian diagonalization approach 
in EIJ. 

It is useful to examine Eqs. (|17H19|) in the limiting (dis- 
sipationless) case for the system by setting j c = Yd = 0. 
In this limit, both coefficients D± and C± become real 
and GLE in Eq. (fTT)) is reduced to a standard nolinear 
Schrodinger equation (NLSE), that possesses exact bright 
and dark soliton solutions depending on the sign of the 



coefficients C± = C^/ and D± = D±' . Here, we suppose 
that coefficients D± > 0. Then, for the case of C± > 0, 
we have a bright soliton solution; while for the case of 
C± < 0, dark solitons can be found 50]. 

In Fig. 3, we show the dependences of this nonlinear 
coefficient, C±\ as a function of the qubit interaction 
strength r\. In Fig. 3, we see that the upper (optical) 



(i) 




FIG. 4. (Color online) Polariton number density n po \ ver- 
sus qubit interaction strength g required for soliton formation 
taken for 7 = 0.1, g = 1. Other parameters are the same as 
in Fig. 3. 



branch gives only a non- negative value of C± ; while the 
lower branch that characterizes matter excitations gives 
only a non-positive value. In other words, we can only 
have a bright soliton solution supported in UB; while 
a dark one in LB. But due to the requirement to con- 
serve the polartion number only bright soliton solutions 
can meet this condition for existing 51]. For a moder- 
ate polariton density number of total excitations, that is 
^poi = 0.01, as shown in the inset of Fig. 3, the value 

of tends to zero for the attractive qubit interaction, 
77 > 0. But to one's surprise, as compared to the posi- 
tive one, a negative qubit interaction 77 < can induce 
a larger nonlinear coefficient, C± , which is enhanced by 
two orders of the magnitude. At the same time, as shown 
in Fig. 2(c), this bright soliton also moves as a "slow 
light" . This giant enhanced nonlinearity is inherent in 
matter-like polaritons and results from the repulsive in- 
teraction among the qubits in this case [l8j |. 

The difference for the nonlinear coefficient C+ between 
positive and negative values of 77 can be understood in 
the following way. For a positive value of 77, the nonlin- 
ear interaction among qubits is attractive. In this region, 
the induced nonlinearity mediated through qubit inter- 
actions is almost zero. Instead, for a negative value of 77, 
the repulsive qubit interaction effectively contributes to a 
great amount of the induced nonlinearity of the photon- 
like UB polariton soliton. Moreover, the properties of 
nonlinear coefficient C± being a function of polariton 



number density n po i are evident from Eq. (|19[) . Effec- 
tively, if one increases polariton number density n po i at 
each cavity, the effect of the ensemble interacting qubits 
results in the enhancement value of a nonlinear coeffi- 
cient. Practically, this effect can be used for creating 
artificial cluster materials based on physically small but 
macroscopical volumes of a qubit-light interaction. Pho- 
tonic crystal defect cavities seem to be an appropriate 
candidate for providing such interaction, cf. [21 1. 
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V. BRIGHT POLARITON SOLITONS UNDER 
PERTURBATIONS 

Now we are going to examine a coupled qubit-light 
system in respect of bright soliton formation; the system 
being non-equilibrium. In this case we are interested in 
soliton solutions of Eq. (fTT|) for optical field amplitudes 
'F in the presence of polariton formation. It is useful to 
bring Eq. (|17[) to a dimensionless form, rewriting it as 
follows 



<9* 1 d 2 ^ 
+ 



dr 2 dx 2 



-iei|*| 2 * + ie 2 



<9 2 * 

dx 2 ' 



(24) 



where 



= X + 1 d and t 



2D\ 
~3T 



■t are new di- 



mensionless spatial and temporal coordinates respec- 



tively; Si 



£2 



|£g| 

2D ( r 



are perturbation coeffi- 



+ 

cients (£1,2 > 0). In Eq. ([24)) we also use the fact that 
the characteristic time scale of a dispersion action and 
the influence of nonlinearity on wave packet spreading 
should be the same [5(|. In particular, the condition 



4z?i 1) In, 



d 2 9 4 



Re 



3n, 



477 



■r 



(25) 




FIG. 5. (Color online) Snapshots of bright soliton profile |$| 
as a function of effective coordinate x e fl = x — f (t) at dif- 
ferent time of propagation r. In the inset temporal behav- 
ior of soliton group velocity v is plotted. The parameters 
are (cf. [H, E|): fl = 2tt x 1.7THz, T = 2tt x 12.1GHz, 
7 C = 2tt x 50GHz, a = 2tt x 0.75THz, 77n po i = 2tt x 24.3GHz, 
d = 400nm, k = 10 4 m _1 , and 5 = 0, respectively. 



should be met for polariton number density n po \ in this 
case. 

Dependences of function of qubit interaction 

strength 77 are presented in Fig. 4. The limiting cases 
for polaritons are determined by the magnitude and sign 
of detuning 5 that includes the dependence on 77. For 
large positive values of 77, low-excitation density limit 
(n po i "C 1) is violated for the chosen qubit-light interac- 
tion parameters. In this case we should take into account 
next terms in expansion of operators J+ ; i characterizing 
a two-level oscillator system, resulting in a quintic non- 
linearity [2l| . 

We examine the GLE in Eq. (p4|) using perturbation 
theory for solitons by consulting the reference In 
particular, we are looking for the solution of Eq. (|24|) in 
the form 



*(r, x) = 2z/sech [2v (x - C(t))] e 1 ^ 



t) 



(26) 



where v, C( T )i an d V 5 ( r ) a; ) are related soliton amplitude, 
position, and phase, respectively. In the absence of per- 
turbation (ei ( 2 = 0), the ansatz used in Eq. (j26|) repre- 
sents an exact solution for Eq. ([24|) with the parameters 
evolving in time as 



VT 

C(r) = T +Co, 

lp (T,x) = ^(x-((T))+5(t), 



Sir) = V 



2v 2 t 



(27) 
(28) 

(29) 



where parameters £0 an d <5o represent the initial soliton 
position and phase respectively. Here, v = v + is soliton 



velocity that we can associate with the velocity of UB 
polaritonic wave-packet by setting 7 = Td = 0. In the 
presence of small perturbation (£1,2 7^ 0), we consider the 
soliton amplitude v and velocity v to be time dependent. 
Applying soliton perturbation theory, we arrive at the set 
of equations for solitons parameters under the adiabatic 
approximation, 



v = -| (2£i + £ 2 ) v 3 - ljS 2 v 2 v, 



v+ 



16 



C=- 
^ 2' 

, v 2 



-E 2 VV 



2v 2 , 



(30) 
(31) 
(32) 

(33) 



where dot denotes derivative in respect of normalized 
time, i.e., d/dr. 

In Fig. 5, we represent a typical UB polariton soli- 
ton profile taken at different time for qubit-light interac- 
tion in a semiconductor GaAs-based quantum well (QW) 
structure, as an example. Based on real parameters of 
GaAs-based quantum wells, the initial shape of this soli- 
ton is characterized by the dashed (red) curve. 

From Fig. 5, it is evident that a soliton amplitude van- 
ishes due to effective nonlinear absorption. At the same 
time the group velocity of soliton is slightly modified 
due to diffusion processes characterized by parameter £2 
shown in Eq. (1301) . The characteristic time of disturbing 
soliton group velocity is tens of picoseconds which seems 
to be reasonable taking into account polariton lifetime in 
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semiconductor structures (the maximal magnitude of di- 
mensionless time r = 12 in the inset of Fig. 5 corresponds 
to the value of 25 picoseconds in real time scale), cf. fl§h 
[20j . Moreover, the temporal soliton displacement £(t) 
can be characterized by Eq. (|2T[) with a good approxi- 
mation. 



VI. CONCLUSIONS 

In this work, we consider the problem of soliton for- 
mation in the array of weakly coupled qubit ensembles 
interacting with the optical field in a tunnel-coupled cav- 
ity array. The effects of cavity field dissipation and qubit 
dephasing have been taken into account. By expanding 
Lipkin-Meshkov-Glick model for the ensemble of inter- 
acting qubits with the cavity photon fields in the array, 
we analytically reveal the existence of mean-field soli- 
ton solutions. In particular, we use multiple-scale en- 
velope function method to obtain an effective complex 
Ginzburg-Landau equation, containing nonlinear absorp- 
tion and diffusion terms. Mean-field soliton solutions 
for this coupled matter-field states are revealed in the 



representation of UB and LB polaritons. With a neg- 
ative value of the qubit interaction, not only a "slow" 
wave-packet in the form of solitons but also an enhanced 
nonlincarity can be achieved in such a qubit-cavity-QED 
array configuration. We have shown that in this case 
an optical wave-packet can be recognized as a macro- 
scopical polariton qubit state. In particular, we have 
demonstrated that bright solitons are formed for UB po- 
lariton wave packets and characterized by slowly varying 
soliton amplitude, momentum, position, and phase. No- 
tably, the group velocity of soliton decreases in time due 
to diffusion processes within a picoseconds domain for 
semiconductor structures. Our analytical and numerical 
results presented here and connected with a cavity-QED 
array provide a platform for the studies of collective prop- 
erties of light with interacting media. 
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